Thermoelectric properties of electrically gated bismuth telluride nanowires 
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We theoretically studied the effect of the perpendicular electric field on the thermoelectric prop- 
erties of the intrinsic, n— type and p— type bismuth telluride nanowires with the growth direction 
[110]. The electronic structure and the wave functions were calculated by solving self-consistently the 
system of the Schrodinger and Poisson equations using the spectral method. The Poisson equation 
was solved in terms of the Newton - Raphson method within the predictor-corrector approach. The 
electron - electron exchange - correlation interactions were taken into account in our analysis. In the 
temperature range from 77 to 500 K, the dependences of the Seebeck coefficient, thermal conductiv- 
ity, electron (hole) concentration, and thermoelectric figure of merit on the nanowire thickness, gate 
voltage, and excess hole (electron) concentration were investigated in the constant relaxation-time 
approximation. The results of our calculations indicate that the external perpendicular electric 
field can increase the Seebeck coefficient of the bismuth telluride nanowires with thicknesses of 7 - 
15 nm by nearly a factor of 2 and enhance ZT by an order of magnitude. At room temperature, 
ZT can reach a value as high as 3.4 under the action of the external perpendicular electric field 
for realistic widths of the nanowires. The obtain results may open up a completely new way for a 
drastic enhancement of the thermoelectric figure of merit in a wide temperature range. 

PACS numbers: 73.63.Nm, 72.20.Pa 
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I. INTRODUCTION 

Nanostructured materials offer a strong potential for 
improving the thermoelectric figure of merit and effi- 
ciency of the thermal-to-electric energy conversion^ 
Bismuth telluride and its solid solutions (Bi2-xSb x Te3, 
Bi2Sb3_j,Se) are among the best known thermoelec- 
tric materials for the modern commercial application. 
These materials possess notable properties such as high 
anisotropic multivalley Fermi surface, small value of the 
thermal conductivity, and optimal value of the carrier 
concentration at room temperature (RT) . Since the elec- 
tron dimensional confinement and the phonon bound- 
ary scattering control the electron and phonon trans- 
port, bismuth telluride based superlattices and quan- 
tum wires have improved thermoelectric properties as 
compared to bulk bismuth telluride material. For ex- 
ample, the thermoelectric figure of merit {ZT) of the 
Bi2Te3/Sb2Te3 superlattices was reported to have the 
value of 2.4 while the maximum ZT for the bulk bismuth 
telluride material is less than 1 at room temperature^ 
Due to stronger spatial confinement the bismuth telluride 
nanowires (NW) are expected to be better than the quan- 
tum well superlattices. Bismuth telluride NWs can be 
fabricated by the electrochemical deposition of the ma- 
terial in the nanopores of anodized alumina membranes, 
by the Taylor-Ulitovsky technique, and by high pressure 
injection of the melt into capillaries Using the tech- 
nique based on translating thin film growth thickness 
control into planar wire arrays, metallic and semicon- 
ductor nanowires with diameters and pitches (center-to- 
center distances) as small as 8 and 16 nm, respectively, 
were prepared^ 



Electric field effect (EFE) is a powerful tool controlling 
the electrical properties of low dimensional structures. 
The latest experimental and theoretical studies showed 
that the EFE, as well as thermal annealing, can signifi- 
cantly improve the thermoelectric properties of Bi, PbSc 
nanowires and PbTc films i 7 ' 8 ' 9 ' While dimensional con- 
finement of electrons leads to the modification of their 
density of states, the Fermi level can be changed due 
to EFE control of the electron (hole) concentration in 
NW. In nanowires, EFE on the thermoelectric proper- 
ties can also be related to the energy spectrum modifi- 
cation and local properties of the electron distribution 
function. In addition, the applied electric field can cause 
quantum confinement of carriers. For these reasons, the 
electric field applied through the side gates camn lead to 
strong modifications the thermoelectric properties of the 
nanowires. The analysis of the temperature dependent 
thermoelectric power provides a reliable way for estimat- 
ing the majority carrier concentration in nanowires even 
when the conventional Hall - effect cannot be used^i 

In this paper, we theoretically studied the intrinsic and 
doped bismuth telluride nanowires grown along [110] di- 
rection at RT. The studied NWs have the square cross- 
sections of 7 x 7 and 15 x 15 nm 2 . The main objective of 
this manuscript is to examine a feasibility of a drastic ZT 
enhancement in nanowires via application of the perpen- 
dicular electrical field. Specifically, we studied how the 
change in the electron density of states caused by both 
the size quantization and electric field influences the ther- 
moelectric properties of NWs. We examined conditions 
at which the size quantization effects compete with the 
electric field effect and when both effects cooperate in 
improving the thermoelectric properties of NWs. The 
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electron mobility in the NWs is assumed to coincide with 
the electron mobility in the bismuth telluride bulk mate- 
rial. In our calculations, we used the values of the lattice 
thermal conductivity, electron (hole) mobility, and car- 
rier effective mass components obtained previously for 
the bismuth telluride NWs with the growth direction of 
[HO]. 1 - The transport parameters were estimated using 
the semi-classical Boltzman approximation. We used the 
constant-relaxation-time approximation to calculate the 
Seebeck coefficient, electrical and thermal conductivities. 
The application of this approach to the calculation of 
the transport properties of bismuth telluride NW ma- 
terial has been justified by us elsewhere^ In this pa- 
per, we also used some of the notations introduced in 
our previous work.— The rest of the paper is organized 
as follows. In Section [IIJ we describe our approach for 
solving the self-consistent system of the Schrodinger and 
Poisson equations. In Section Hm the dependence of the 
thermoelectric parameters of bismuth telluride NWs on 
lateral gate voltage is studied. In Section IIV1 we con- 
sider the local variation of the thermoelectric parameters 
inside NWs. In SectionfV] we discuss the temperature de- 
pendence of the thermoelectric parameters. Finally, our 
conclusions are given in Section IVII In the Appendix, a 
detailed description of the spectral method, based on the 
trigonometric functions, and the spectral element method 
(SEM), based on the Gauss-Lobatto-Lelendre polynomi- 
als, and the linearized Poisson equation (LPE) approach 
are presented. 



II. THEORY 



A. Self-consistent system of the Schrodinger and 
Poisson equations 



We considered a bismuth telluride nanowire connected 
to the source and drain and separated from the lateral 
metallic gate by the dielectric SiC>2 layer with the thick- 
ness e?si0 2 = 300 nm. The y axis of the system of coor- 
dinates is directed along the nanowire axis. The elec- 
tronic band structure of the indirect band gap bismuth 
telluride NW material is given by a parabolic band ap- 
proximation in the Drable-Wolf six- valley model.— In the 
rhombohedral coordinate system, the main axis of two 
equivalent electron (hole) ellipsoids is along the direction 
[110]. For the other 4 electron (hole) equivalent ellip- 
soids, the main axis is at the angle 7r/3 with respect to 
the direction [110]. The applied electric field is assumed 
to be parallel to the z axis. A similar device based on the 
Bi nanowire with a thickness of 28 nm has been recently 
studied experimentally^ To estimate the nanowire trans- 
port properties, we found a solution of the self-consistent 
system of the Schrodinger and Poisson equations taking 
into account the difference between the effective masses 
for the non-equivalent ellipsoids 
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+ [V c (x, z) + V xc (x, z)+E c - e<p(x, z)} z) = E^(x, z), 



(1) 



ee Aip(x, z) = p(x,z), 



(2) 



where the local variables x and z range in the inter- 
vals < x < a x and < z < a z , E c is the bottom of the 
conduction band, (p is the electrostatic potential, and 
V c is the confinement potential, which corresponds to 
an infinitely high quantum well. The charge density 



is given by the expression p = epiu — enuj, where nm 
(jPxd) is the electron (hole) concentration. V xc {x,z) is 
the electron-electron exchange-correlation energy given 
by the following expressio n 13 i 14 



V xc (x, z) 



1 + 0.03683r s ln 1 



21 

r s (x,z). 



j s (x, z) na ) Ry 



(3) 



where a — (47T/9) 1 / 3 , a* = 4ire eh 2 / 'e 2 m* is the Bohr ra- 
dius, r s (x, z) = (3/47rn (x, z)) 1 / 3 /a*, R* = e 2 /87re ea* 
is the effective Rydberg constant. The dielectric con- 
stant of bismuth telluride is esiTe = 90 (RefllU8lll5l). In 



the exchange-correlation potential, the anisotropic effec- 
tive mass m* can be chosen as either a constant av- 
erage mass m* c = (m x m y 'm z } 1 l z or an optically aver- 
aged mass m* c = 3/(l/m x + l/m y + l/m z ) . For both 
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choices of the effective mass, the results of energy cal- 
culation are similar^ For the sake of simplicity, in our 
calculations we use the constant average mass. For the 
NW growth direction [110], the electron constant aver- 
age mass is m% c = 0.0821 toq. The hole average mass is 
m h xc = 0.105 mo- Due to the small values of the electron 
(hole) effective masses and a large value of the dielectric 
constant, the Bohr radius is a* = 57 nm and a* = 45 
nm for electrons and holes, respectively. It is greater by 
3 orders of magnitude than the Bohr radius of hydrogen 
a* = 0.053 nm. For the electrons (holes), the Rydberg 
constant is R* = 1.37 x 10~ 4 eV (R* y = 1.77 x 10" 4 eV) 
which is nearly 2 orders of magnitude less than the Ry- 
dberg constant for hydrogen R* — 1.1 x 10 -2 eV. Since 
the dielectric constant of bismuth tclluride material is 
high, the Rydberg constant has a rather small value of 
about 15 meV compared to the band gap of 480 meV 
for the NW with a thickness of 7 nm and the excess 
hole concentration p ex = 1 x 10 19 cm~ 3 . The influence 
of the electron-electron exchange-correlation processes 
on the higher-temperature performance of the wires is 
largely irrelevant, as the values of V xc do not exceed 20 
meV, which corresponds to about 230 K. For the in- 
trinsic semiconductor nanowires with a small electron 
(hole) concentration, we did not take into consideration 
the local exchange-correlation energy, since the exchange- 
correlation effects can be neglected in semiconductor het- 
erostructures with the doping concentration of about 
10 18 cm -3 (Ref. [l7T ). After switching the sign in the 
second derivatives in Eq. (TT|), we obtain the Schrodinger 
equation for the holes. 

Integrating both sides of Poisson's equation, we obtain 
the charge neutrality equation for the gated NW in the 



n 1D (E F ) = pid(E f ) + a x a. 



(4) 



Here, the parameter a x denotes the NW thickness along 
the x axis. The surface electric charge — ea , induced 
by the gate electric field at the Si02/Bi 2 Te3 interface, 
is associated with a Fermi-level-pinning boundary con- 
dition. The condition ^ means that the number of 
electrons and holes that populate the subbands is equal 
to the number of charge carriers, which occupy the sur- 
face states. The surface states do not contribute to the 
electron transport, because of the strong surface rough- 
ness scattering. The applied electric field induces a sur- 
face charge on the SiC>2/Bi2Te3 interface. The induced 
surface charge doping concentration a can be estimated 
in the gate — SiC>2 — Bi2Te3 capacitance structure as a 
function of the gate voltage V g (Ref. [pj) 



ed s 



Va 



i0 2 



(5) 



where esi0 2 = 3-9 is the permittivity of the SiC>2 
layer. The effective doping induced by the ap- 
plied electric field can be estimated from the expres- 
sion Nd = V g eoesio 2 / e dsi0 2 a NW , where clnw is the NW 



thickness. In the above expression, the gate volt- 
age V g is followed by the factor of 4.79 x 10 16 cm~ 3 /V^ 
(1.02 x 10 17 cm~ 3 /V) for the NW thickness 15 nm (7 
nm). Hence, the NW electron (hole) concentration can 
be adjusted to the desirable value by tuning the gate volt- 
age. The same effect can be achieved by doping NWs. 
However, it is accompanied by an undesirable increase of 
the electron-impurity scattering. Since the electric field 
decreases across the dielectric layer, the effective electro- 
static potential V* at the Si0 2 /Bi 2 Te 3 interface is given 
by the expression V* = V g e Sl o 2 aNW 1 1 tBi 2 Te 3 d Sl o 2 - The 
effective gate voltage V* at the Si02/Bi 2 Te3 interface is 
less than the gate voltage V g by a factor of 460 (10 3 ) 
for the NW thickness of 15 nm (7 nm). The boundary 
conditions for Poisson's equation are 



tp(z = 0) = 0, 



ea. 



(6) 



(7) 



To solve the Schrodinger equation, we use the spectral 
element method (SEM),! 9 .! 2 ^ 2 ! For the intrinsic NWs, 
we solve Poisson's equation after its linearization. 22 In 
this case, the Schrodinger equation can be recasted 
into the Mathieu equation. 23 i 24 For the doped NWs, we 
solve the nonlinear Poisson equation by means of the 
Newton-Raphson method using the predictor-corrector 
approachj 25 i 26 i 27 In the next subsection, we describe our 
method of calculation of the thermoelectric transport co- 
efficients taking into account the local and gate voltage 
dependences of both the electron wave functions and the 
energy spectrum. 



B. Thermoelectric transport coefficients 

To calculate thermoelectric parameters, we use the 
constant relaxation time (r = To) approximation. In this 
approximation, the electron (hole) mobility [i is assumed 
to be constant and given as [i = er /m*. The electri- 
cal conductivity a , Seebeck coefficient S, and thermal 
conductivity K e (h) f° r a NW are defined a o 12 ' 28 



- r.(°) 



a = L 



S = - 



1 lw 

eT m ' 



(8) 



(9) 



(10) 



The lattice thermal conductivity measured perpendicu- 
lar to the trigonal axis, kl , is equal to 1.45 W/mK for 
bulk bismuth telluride at RT ; 29 ' 30 It is inversely propor- 
tional to temperature, n L = 430/T, when T > 50 KM 
Both the theoretical estimations and measurements of 
the NW lattice thermal conductivity indicate that it is 
reduced by an order of magnitude compared to its bulk 
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counterpartj 32 ' 33 We take into consideration this fact in 
the calculation of ZT. The gate electric field does not in- 
fluence the lattice thermal conductivity. The transport 
matrix elements are defined below in Eqs. p2[) and 
(fl~3|) . where the electron (hole) part of the thermal con- 
ductivity is given by 



matrix element includes both the electron and hole 



parts, = Li a) + L^ } . In the LPE approach, the 
form of the transport matrix elements coincides with 
that given for the bismuth telluride NWs in the absence 
of the gate voltage^ 



-(«) 



e(h) 



The 



last 



1 

term 



L (2)_ 



(ii) 



Eq. (nnD, 

attributed to 



K eh = Ta e a h (S h - S e ) 2 /(a e + <j h ), is 
the thermal energy carried by the electron-hole pairs 
generated at the heated end of the sample and moving 
to the cold end where they annihilate. The generalized 
transport matrix element L' a ) represents the sum 
over all six valleys of the bismuth telluride material, 

l(q) = Ei L i a) - For the bipolar nanowire, the transport 



L? 



E 



dk„ 



vf{k y ) [E n j{k y ) - E F ] a dk y . 

(12) 



In this case, the electric field effect on the thermoelec- 
tric parameters is due to the modification of the energy 
spectrum caused by the electrostatic potential ip(x, z). In 
the framework of SEM, the local variation of the trans- 
port parameters is defined by the electron (hole) wave 
function^ In SEM approximation, the transport matrix 
elements for each valley have the form 



Lf(x,z) 



Tra x a, 



- 1 

n.L 



(13) 



r 



In this case, the electron concentration for each valley 
can be written as 

n 1D (x, Z )=Nl° ®-l/2(Vn,l)\Vn,l( X ,z)\ 2 - 
n=l 1=1 

The thermoelectric parameters are effectively macro- 
scopic parameters. Hence, under an applied lateral elec- 
tric field, the figure of merit ZT and the Seebeck coeffi- 
cient S are defined by averaging as 8 



{a{z)S{z)Y 
(v(z))(k(z)) ' 



(15) 



(16) 



In the next section, we present the results of the cal- 
culation of thermoelectric parameters in dependence on 
gate voltage in terms of the Spectral Element Method 
and the Linearized Poisson's Equation approach at room 
temperature. 



III. DEPENDENCE OF THE 
THERMOELECTRIC PARAMETERS ON THE 
GATE VOLTAGE 

For the intrinsic NWs, our calculations show a good 
agreement between the results obtained in the LPE and 



SEM approximations. Figure [T] presents the depen- 
dence of the Debye screening length Ljj for the intrinsic, 
n— type, and p— type bismuth telluride NWs with thick- 
nesses of 7 and 15 nm on the gate voltage V g at room 
temperature. The Debye screening length decreases with 
decreasing NW thickness. The maximum value of the De- 
bye screening length is achieved at V g = and V g = 50 V 
(-50 V) for the intrinsic and p— type (n— type) NWs, cor- 
respondingly, with a thickness of 7 nm when the carrier 
concentration is minimal. There are no peaks in the gate 
voltage dependence of the type (n— type) NWs with 
a thickness of 15 nm, because their type of conductivity 
is not changed in the given range of gate voltage. The 
gate voltage dependence of the Debye screening length 
for the 15-nm-thick NW is much less than that for the 
7-nm-thick NW. For the 15-nm-thick intrinsic NW, Ljj 
decreases nearly twice. It is less than the NW thickness 
in the whole range of the gate voltage. Lo decreases 
by 9 times for the intrinsic NW with a thickness of 7 
nm. It is almost equal to the NW thickness at V g = ±30 
V. The ratio of the maximal value of Ljj to its mini- 
mal value is 6 and 1.3, respectively, for the doped NWs 
(Pex = 5 x 10 18 cm -3 ) with thicknesses of 7 and 15 nm. 
Therefore, the EFE allows controlling the value of the 
Debye screening length for both the intrinsic and doped 
bismuth telluride NWs. 

Figure [2] depicts the dependence of the band diagram 
of the intrinsic and p— type bismuth telluride NWs with 
thicknesses of 7 and 15 nm on applied gate voltage with 
respect to the Fermi energy (Ep = 0). The calculations 
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FIG. 1: (Color online) Dependence of the Debye screening length on gate voltage (a) for the intrinsic bismuth telluride NWs 
with the thicknesses a z = 7 (solid line) and 15 nm (dashed line), (b) for the doped bismuth telluride NWs with thicknesses of 
7 (solid line) and 15 nm (dashed line) at the excess hole concentration p ex — 5 x 10 18 cm -3 . The dashed-dotted line (dotted 
line) corresponds to the NW with a thickness of 7 nm and excess hole concentration p ex = 1 x 10 19 cm -3 (excess electron 
concentration n ex = 5 x 10 18 cm -3 ). The data obtained by using the LPE approach are depicted by the filled squares and 
triangles. 




FIG. 2: Dependence of the electron and hole spectrum on gate voltage (a) for the intrinsic and (b) for the p— type bismuth 
telluride NWs with the thicknesses a z = 7 nm (solid line) and 15 nm (dashed line), when the excess hole concentration is 
Pex — 5 x 10 18 cm -3 at room temperature. 



TABLE I: Modification of the electronic structure for the intrinsic bismuth telluride NWs under influence of the gate voltage 
at room temperature. 

Thickness, nm V g , V AEf 2 , meV A£f 3 , meV AE? 2 , meV AE% 3 , meV AE g , meV 

7 109 33 73 5 322 

7 70 102 38 66 10 340 

15 23 7 16 1 171 

15 70 35 10 25 16 160 
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are based on the SEM. Due to the greater size quan- 
tization effect, both the separation between the energy 
subbands and the band gap is greater for the NW with a 
thickness of 7 nm than those for the NW with a thickness 
of 15 nm. At the positive gate voltage, the bottom of the 
conduction band for the NW thickness of 7 nm coincides 
with that for the NW thickness of 15 nm. At the negative 
gate voltage, the top of the valence band for both NW 
thicknesses is the same. This approaching of the sub- 
bands of the NWs with different thicknesses reflects the 
electric field doping effect, which changes the NW polar- 
ity, when the Fermi level approaches the top (bottom) 
of the valence (conduction) subband. Since, for the NW 
thickness of 7 nm, the Debye screening length is greater 
than the NW thickness at the small values of the gate 
voltage \V g \ < 25 V (see Fig. [T]), the dependence of both 
the electron and hole energy on gate voltage is stronger 
than that for the NW thickness of 15 nm. For the 7-nm- 
thick (15-nm-thick) NW, the band gap increases nearly 
by 5 per cent from 323 meV (172 meV) to 340 meV (183 
meV), when the gate voltage ranges in the interval from 
to 70 V (25 V). For the 15-nm-thick NW, the behavior 
of the band gap dependence is not monotonic because, 
at a large value of the gate voltage, the applied elec- 
tric field influences the size quantization energy levels of 
electrons and holes in different manner due to the nonlin- 
ear electrostatic potential. For the 15-nm-thick NW, the 
electrostatic potential f(V g ) behaves as a parabolic func- 
tion of gate voltage, f(V g ) — aV g 2 + bV g + c, where a, b, 
and c are some constants. It behaves as a linear function 
of gate voltage, f(V g ) — aV g + b, for the 7-nm-thick NW. 
This different behavior is due to the fact that the Debye 
screening length for the 7-nm-thick NW is about equal to 
the NW thickness, while it is about twice less than the 
thickness of the 15-nm-thick NW (see Fig. [TJ. Table U 
shows that, under the influence of the electric field, the 
electronic structures for the intrinsic NWs with thick- 
nesses of 7 and 15 nm are modified in different manner. 
For the intrinsic 7-nm-thick (15-nm-thick) NW, the 1st, 
the 2nd, and the 3rd electron subbands move up by 149 
(74), 156 (62) and 151 (59) meV, correspondingly, when 
the gate voltage varies from to 70 V. For the 7-nm-thick 
NW, the spacing between the 1st and the 2nd electron 
(hole) subbands A_Ef 2 (AE^ l 2 ) decreases by 7 meV, while, 
for the 15-nm-thick NW it increases by 8 meV (9 meV). 
Also, for the n— and p— type NWs, the modification of 
their electronic structure depends on the NW thickness. 

Figure [3] shows the linear dependence of the electron 
and hole concentrations in the intrinsic and p— type bis- 
muth telluride NWs on the gate voltage at room tempera- 
ture. At the zero gate voltage, the charge carrier concen- 
trations in the intrinsic bismuth telluride NWs with the 
thicknesses of 7 and 15 nm are n = p = 3.11 x 10 16 cm" 3 
and n = p = 3.03 x 10 17 cm" 3 , respectively. Under the 
applied lateral electric field, the carrier concentration in- 
creases by two orders of magnitude for the 7-nm-thick 
intrinsic NW and by one order of magnitude for the 15- 



nm-thick intrinsic NW. Since the hole effective mass is 
greater than the electron effective mass, the hole concen- 
tration increases faster with the gate voltage compared 
to the electron concentration. For this reason, for the 
intrinsic NW with a thickness of 15 nm, the dependence 
of the electron concentration on the gate voltage is less 
pronounced compared to that for the holes. At the pos- 
itive (negative) gate voltage, the electrons (holes) domi- 
nate. Therefore, it is possible to control the type of the 
nanowire conductivity by applying gate voltage. This 
fact was experimentally confirmed for Bi nanowires^ At 
the excess hole concentration p ex = 5 x 10 18 cm" 3 , the 
hole concentration can increase nearly twice and three 
times for the type NWs with thicknesses of 7 and 15 
nm, respectively. The charge carrier concentration di- 
minishes, while the gate voltage dependence of the car- 
rier concentration rises up with a decrease in the NW 
thickness. This effect arises from the increased splitting 
between the energy subbands as a result of a large con- 
finement effect. Thus, the number of subbands involved 
in the transport decreases with a decrease in the nanowire 
thickness. For the n— type NW, the dependence of the 
carrier concentration on the gate voltage is reversed com- 
pared to that for the p— type NW. For the 7-nm-thick 
NWs with different excess hole concentrations, the de- 
pendences of the carrier concentration on the gate volt- 
age are the same. The effective doping Nd induced by the 
applied electric field is proportional to the gate voltage 
and inverse proportional to the NW thickness (see dis- 
cussions below Eq. [5]). By proper adjusting of the gate 
voltage, the effective electric doping can compensate the 
excess holes (electrons). Hence, the point of neutrality 
(Vg) for the p(n)-type NW, at which the NW is getting 
to be intrinsic, depends on both the excess hole (elec- 
tron) concentration and NW thickness. Therefore, the 
point of neutrality (V g ~ 50 V) for the 7-nm-thick p-type 
NW is less than that (V g w 100 V) for the 15-nm-thick 
NW with the same excess hole concentration. The gate 
voltage dependence of the carrier concentration rises up 
with a decrease in the NW thickness due to the strong 
size quantization effect. Hence, the carrier concentration 
corresponding to the 15-nm-thick NW is greater (less) 
than that for the NW with a thickness of 7 nm at the 
V g = +50 V (V g = -50 V). 

In our calculations, we assumed that the carrier mobil- 
ity is constant. Since the electrical conductivity is pro- 
portional to the carrier concentration, the dependence of 
the electrical conductivity on the gate voltage Vg is sim- 
ilar with that for the carrier concentration depicted in 

Fig. m 

Taking into account the gate voltage dependence of 
the carrier concentration, we calculated the exchange- 
correlation energy as a function of gate voltage and NW 
thickness. For the p— type NW with a thickness of 7 
nm and excess hole concentration p ex = 5x 10 18 cm" 3 
(p ex = 1 x 10 19 cm~ 3 ), the exchange-correlation energy 
changes from the value V xc = — 10 meV (-13 meV) at the 
gate voltage V g = to V xc — — 15 meV (-17 meV) at the 
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FIG. 3: (Color online) Gate voltage dependence (a) of the electron (n) and hole (p) concentrations for the intrinsic bismuth 
telluride NWs with the thicknesses a z — 7 nm (solid line) and 15 nm (dashed line), (b) of the electron and hole concentrations 
for the p— type NWs with thicknesses of 7 nm (solid line for holes and dashed-dotted line for electrons) and 15 nm (dashed 
line for holes) when p ex — 5 x 10 18 cm -3 , at room temperature. The dotted line corresponds to the gate dependence of the 
hole concentration for the NW thickness of 7 nm when p ex = 1 x 10 19 cm -3 . The data obtained in the LPE approximation are 
depicted by the filled squares and triangles. 



gate voltage V g — —85 V. At the excess hole concentra- 
tion p ex — 5 x 10 18 cm" 3 , the exchange-correlation en- 
ergy does not change with increasing gate voltage, while, 
at the excess hole concentration p ex = 1 x 10 19 cm" 3 , it 
achieves the value V xc — —6 meV at the gate voltage 
V g = 85 V. For the p— type NW with a thickness of 15 nm 
and the excess hole concentration p ex = 5 x 10 18 cm" 3 , 
the exchange-correlation energy varies from the value 
V xc = — 16 meV to the values V xc = —7 meV and -20 
meV at the gate voltages V g = 85 V and V g = —85 
V, correspondingly. For the n— type NW with a 
thickness of 7 nm and the excess electron concentra- 
tion p ex — 5 x 10 18 cm" 3 , the exchange-correlation en- 
ergy varies from the value V xc = — 10 meV to the value 
Vxc = — 15 meV at the gate voltage V g = 85 V, while V xc 
does not change with decreasing gate voltage. 

Figure @] shows the dependence of the nanowire Seebeck 
coefficient (S) on gate voltage at room temperature. At 
the zero gate voltage, the Seebeck coefficient is -167 and 
-252 V/K for the intrinsic NWs with thicknesses of 15 
and 7 nm, correspondingly. For the intrinsic NW with 
a thickness of 7 nm, the maximum absolute value of the 
Seebeck coefficient approaches the value of 478 V/K (494 
V/K) at the gate voltage of -3 V (2 V). For the 15-nm- 
thick intrinsic NW, the maximum absolute value of 333 
V/K (344 V/K) is achieved at the greater gate voltage of 
-29 V (15 V). Since the electron mobility, 1200 cm 2 /Vs, 
in the binary-bisectrix plane of the bismuth telluride ma- 
terial is twice as large as the hole mobility, 510cm 2 /Vs 
(Ref. fl2h . the maximum absolute values of the Seebeck 
coefficient corresponding to the opposite gate polarities 



are achieved at different values of V g . For the intrinsic 
NWs, the EFE causes the Seebeck coefficient to increase 
as much as twice. The confinement effect on the elec- 
tron density of states leads to an increase in the NW 
Seebeck coefficient, while the increase in the carrier con- 
centration causes a decrease in S (Ref. IT2T ). The larger 
electron density of states D(E) for NWs means a small 
Fermi-Dirac distribution function f{E c — Ep) and thus 
a larger difference (E c — Ep). As a result, the Seebeck 
coefficient for nanowires is higher than that of their bulk 
counterparts, whereas the electrical conductivity is the 
same. Therefore, the absolute value of the Seebeck coef- 
ficient strongly increases at the low gate voltage due to 
the increase in the confinement effect. If the gate voltage 
continues to increase, IS"! decreases because the carrier 
concentration increases with the gate voltage. In the gate 
voltage interval from -90 to 90 V, the p— type NW with 
a thickness of 15 nm (7 nm) and the excess hole concen- 
tration p ex = 5 x 10 18 cm" 3 (p ex = 1 x 10 19 cm" 3 ) does 
not change its type of conductivity (see Fig.[3]Jb)); hence, 
the dependence of the Seebeck coefficient of such NW on 
the gate voltage is monotonic. It is increased by 60 per 
cent (by 1.3 times) with gate voltage. For the 7-nm-thick 
p— type (n— type) NW with the excess hole (electron) 
concentration 5 x 10 18 cm~ 3 , the absolute value of the 
Seebeck coefficient \S\ varies from the value of 187 V/K 
(206 V/K) at the zero gate voltage to the value of 523 
V/K (507 V/K) at V g = -50 V (50 V) when the NW type 
of conductivity changes (see Fig. [3](b) ) . In this case, the 
EFE increases |5| by 2.8 (2.5) times. Thus, both the con- 
finement effect and the electric field effect considerably 
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FIG. 4: (Color online) Gate voltage dependence of the Seebeck coefficient (S) (a) for the intrinsic and (b) for the p— type 
bismuth telluride NWs with the thicknesses a z = 7 nm (solid line) and a z = 15 nm (dashed line) when p ex — 5 x 10 18 cm -3 
at room temperature. The dashed-dotted line (dotted line) corresponds to the S for the n— type (p— type) bismuth telluride 
NW with a thickness of 7 nm when n ex = 5 x 10 18 cm -3 (p ex = 1 x 10 19 cm -3 ). The data obtained in the LPE approach are 
depicted by the filled squares and triangles. 



improve the Seebeck coefficient of NW. 

Figure [5] depicts the gate voltage dependence of the 
thermoelectric conductivity for the NWs with thicknesses 
of 7 and 15 nm. For the intrinsic bismuth telluride NW 
with a thickness of 15 nm, a small disagreement between 
the LPE and SEM data appears owing to the great num- 
ber of the subbands involved in the transport and the 
rather complex mathematical expression corresponding 
to the electron part of the thermal conductivity. The de- 
pendence of the intrinsic NW thermal conductivity on the 
gate voltage is not monotonic. The confinement effect on 
the electron density of states leads to an increase in both 
the NW Seebeck coefficient and the thermal conductiv- 
ity, while the increase in the carrier concentration causes 
their decreased At the zero gate voltage, the thermal 
conductivity is 0.248 W/Km (0.702 W/Km) for the NW 
with a thickness of 7 nm (15 nm). The minimum value 
of the NW thermal conductivity of 0.18 W/Km (0.39 
W/Km) is achieved at the gate voltage of 2 V (-50 V). 
For the intrinsic NW with thicknesses of 7 nm and 15 
nm, the EFE decreases the thermal conductivity by 30 
and 45 per cent, respectively. The interval between the 
extreme points of the gate voltage for the 7-nm-thick NW 
is less than that for the 15-nm-thick NW because of the 
strong dependence of the carrier concentration on gate 
voltage for the NW with small thickness (see Fig. [3]) . For 
the 7-nm-thick NW of the p— type (n— type) at the ex- 
cess hole (electron) concentration of p ex = 5 x 10 18 cm -3 
(nex = 5 x 10 18 cm -3 ), the minimal thermal conductiv- 
ity of 0.152 W/Km (0.147 W/Km) is achieved at the gate 
voltage of -50 V (50 V) when the type of the NW elec- 



trical conductivity changes into the opposite one. There- 
fore, the EFE halves the thermal conductivity of the 7- 
nm-thick (15-nm-thick) p— type NW when applying gate 
voltage. The difference between the extreme values of the 
p— and n— type NW thermal conductivities is due to the 
difference between the electron and hole effective masses 
as well as the difference in their mobility. For the 7-nm- 
thick (15-nm-thick) p— type NW at p ex = 5 x 10 18 cm -3 , 
the gate voltage dependence of the thermal conductivity 
is monotonic, because the type of NW electrical conduc- 
tivity does not change in the given gate voltage interval. 

Figure [5] presents the gate dependence of the thermo- 
electric figure of merit (ZT) at room temperature. At the 
zero gate voltage, for the intrinsic NWs with a thickness 
of 7 nm (15 nm), the figure of merit is 0.065 (0.099). At 
the gate voltage of 15 V (42 V), the maximum value of 
the figure of merit is as high as 3.4 (2.4). For the opposite 
gate polarity, when the holes dominate, the extreme value 
of the figure of merit is lower than that corresponding 
to the positive gate voltage because of the difference be- 
tween the hole mobility and the electron mobility as well 
as due to the difference in the effective masses of electrons 
and holes. For the 7-nm-thick n— and p— type NWs, the 
maximum value ZT = 3.1 is less than that value for the 
intrinsic NW, ZT = 3.4. For the 7-nm-thick NW, both 
the figure of merit and the Seebeck coefficient have two 
extreme values due to the change of the NW conductiv- 
ity type. When the electrons dominate in the 7-nm-thick 
P~ type NW, the figure of merit of 3.1 at V g = 60 V is 
greater than its value of 1.97 at V g = 25 V by 60 per 
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FIG. 5: (Color online) Dependence of the thermal conductivity on gate voltage (a) for the intrinsic and (b) for the p— type 
bismuth telluride NWs with thicknesses of 7 nm (solid line) and 15 nm (dashed line) when p ex = 5 x 10 18 cm -3 at room 
temperature. The dashed-dotted line (dotted line) corresponds to the thermal conductivity of the n— type (p— type) bismuth 
telluride NW with a thickness of 7 nm when n ex = 5 x 10 18 cm" 3 (p ex — 1 x 10 19 cm" 3 ). The data obtained in the LPE 
approximation are depicted by the filled squares and triangles. 



cent (see Fig. EJb)). At the room temperature, for the 
7-nm-thick p— type NW, the maximum value of the fig- 
ure of merit 3.1 is greater than the corresponding value 
ZT = 1.65 at the zero gate voltage by 90 per cent. Since, 
for the 15-nm-thick NW at p ex = 5x 10 18 cm -3 and the 
7-nm-thick NW at p ex = 1 x 10 19 cm -3 , the type of the 
electrical conductivity does not change within the whole 
interval of the gate voltage, their figure of merit gradually 
increases from 1.19 and 0.93 at the zero gate voltage to 
1.48 and 1.97 at V g = 45 V and V g = 75 V, respectively. 
Therefore, the applied lateral electric field improves the 
NW thermoelectric properties considerably. 



IV. LOCAL VARIATION OF THE 
THERMOELECTRIC PARAMETERS 

In this section, the local variation of thermoelectric 
parameters corresponding to the maximum value of the 
figure of merit at room temperature is discussed. The 
calculations are based on the SEM. For the gated bis- 
muth telluride NWs, the local (maximal and minimal) 
values of thermoelectric parameters as well as their val- 
ues averaged over the NW cross-sections are listed in Ta- 
ble HU The electron (hole) concentration tends to zero at 
the NW boundary, because we impose the zero boundary 
condition on the electron (hole) wave function here. For 
the 7-nm-thick (15-nm-thick) intrinsic NW, the averaged 
value of the electron concentration is 1.5 x 10 18 cm -3 
(1.9 x 10 18 cm" 3 ) at the gate voltage V g = 15 V [V g = 42 
V). The electron distribution in the NW cross section 



is more gradual for the small NW thickness because of 
the strong size quantization effect and the small num- 
ber of subbands participating in the transport for such 
NWs. As a result of the strong size quantization, for the 
less NW thickness, the greater gate voltage should be 
applied to modify the electron dimensional probability 
density function. For the 7-nm-thick (15-nm-thick) in- 
trinsic NW, the Debye screening length 9.5 nm (8.2 nm) 
is less (greater) than the NW thickness at a gate volt- 
age V g = 15 V (Vg = 42 V) (see Fig.CQa)). It also causes 
the uniform (non-uniform) charge carrier distribution in 
the 7-nm thick (15-nm-thick) NW cross section. For the 
n— and p— type NWs with a thickness of 7 nm, their 
optimal profiles of the electron distribution in the NW 
cross section are similar at the gate voltages V g = — 35 V 
and V g — 60 V, when the figure of merit is maximal (see 
Fig.[5{b)). The corresponding Debye screening length, 10 
nm (11 nm), is greater than the thickness of the n— type 
(p-type) NW (see Fig. QJb)). For the n-type NW, 
the optimal profile of the charge carrier distribution is 
obtained at less gate voltage compared to that for the 
p— type NW. Hence, the gate voltage enables charge car- 
rier concentration for the NWs to be manipulated. The 
local variation of the exchange-correlation energy is sim- 
ilar to that of the charge carrier concentration. 

For the intrinsic NW with a thickness of 7 nm (15 nm) 
at V g — 15 V (V g — 42 V), the averaged value of the See- 
beck coefficient is -324 V/K (-306 V/K). For the intrinsic 
15-nm-thick NW, the minimal value of \S\, 257 V/K, is 
achieved at the local point corresponding to the extreme 
value of the NW electron concentration (see Fig. EJa)). 
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FIG. 6: (Color online) The gate voltage dependence of the figure of merit (ZT) (a) of the intrinsic and (b) of the p— type 
bismuth telluride NWs with the thicknesses a z — 7 nm (solid line) and a z = 15 nm (dashed line) when p ex = 5 x 10 18 cm -3 at 
room temperature. The dashed-dotted line (dotted line) corresponds to the ZT of the n— type (p— type) bismuth telluride NW 
with a thickness of 7 nm when n ex — 5 x 10 18 cm -3 (p ex = 1 x 10 19 cm -3 ). The data obtained in the LPE approximation are 
depicted by the filled squares and triangles. 



TABLE II: Thermoelectric parameters of the gated bismuth telluride NWs at room temperature. 



Intrinsic 



p— type 



n— type 



Thickness, nm 

v g , v 

Pex, n ex , 10 18 cm" 3 
Ld, nm 

Umax, 10 18 Cm" 3 

< n >, 10 18 cm" 3 
n mi „ , IlH cm -3 



Pmax, 10 18 Cm" 3 



10 1 
I 1 

< p >, 10 18 cm -3 

Pmm, 10 18 Cm" 3 
Smax, 10" 6 V/K 

< s >, icn 6 v/k 

S mln , 10" 6 V/K 

W/mK 
< k >, W/mK 
W/mK 
< ZT > 



7 
15 

9.5 
2.2 
1.5 
0.53 
8.2 x 10" 3 
5.6 x 1(T 4 

340 
324 
320 
0.313 
0.273 
0.204 
3.4 



15 
42 

8.5 
5.6 
1.9 

0.64 
0.025 


381 
306 
257 
0.704 
0.426 
0.323 
2.4 



7 

60 
5.0 
11 
1.9 
1.2 
0.21 



322 
319 
318 
0.272 
0.223 
0.160 
3.1 



15 
45 
5.0 
7.6 



7.3 
2.8 
0.31 
374 
242 
213 
0.453 
0.278 
0.172 
1.5 



7 

75 
10 



4.4 
2.3 
0.25 
296 
276 
274 
0.280 
0.217 
0.157 
2.0 



7 
-35 
5.0 

10 
2.1 
1.4 
0.39 



304 
301 
300 
0.286 
0.240 
0.173 
3.1 



When the electron and hole local concentrations are the 
same, the corresponding absolute value of the local See- 
beck coefficient, 381 V/K, is maximal. For the 7-nm- 
thick intrinsic and doped NWs with the excess electron 
(hole) concentration 5 x 10 18 cm -3 , the local variation 
of the Seebeck coefficient is weak, which is coherent with 
the local variation of the NW carrier concentration. The 



Seebeck coefficient is negative for both the p— type and 
the n— type bismuth telluride NWs with a thickness of 7 
nm when p ex — 5 x 10 18 cm -3 and n ex = 5 x 10 18 cm" 3 , 
respectively, at the applied gate voltages 60 and -35 V. 
Whereas, S is positive for the p— type NWs with thick- 
nesses of 7 and 15 nm when p ex = 1 x 10 19 cm" 3 and 
Pex — 5 x 10 18 cm" 3 , correspondingly, at V g = 75 and 45 
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V. For the doped NWs, S achieves its maximal value 
when the carrier concentration tends to its minimal value. 
Hence, the applied lateral electric field allows controlling 
the profile of the Seebeck coefficient in the cross section 
of the 15-m-thick NW. 

For the 7-nm-thick NW, the local variation of the 
thermal conductivity is more gradual compared to that 
for the 15-nm-thick NW because of the same reasons 
described above that justify the gradual local varia- 
tion of the charge carrier distribution. For the 7- 
nm-thick n— and p— type NWs, at the excess elec- 
tron and hole concentrations n ex = 5x 10 18 cm" 3 and 
p ex = 1 x 10 19 cm" 3 , respectively, at the gate voltages 
V g = —35 V and V g — 75 V, the local variations of their 
thermal conductivities are the same. As a result, the av- 
eraged values of their thermal conductivities, 0.24 and 
0.22 W/mK, are very close. The increase in the local 
thermal conductivity is related to the increase in the 
charge carrier concentration at the given point of the 
NW cross section. The thermal conductivity is less for 
the NW thickness of 7 nm than that for the NW with a 
thickness of 15 nm due to the enhancement of the dimen- 
sional confinement with a decrease in the NW thickness. 
Therefore, the electric field effect permits controlling pro- 
file of the thermal conductivity in the NW cross section. 



V. TEMPERATURE DEPENDENCE OF THE 
THERMOELECTRIC PARAMETERS 

In this section, we consider the temperature depen- 
dence of the thermoelectric parameters of the 7- and 15- 
nm-thick bismuth telluride NWs at the gate voltages cor- 
responding to a maximum value of the figure of merit at 
room temperature. 

Figure [7] presents the temperature dependence of the 
first three electron and hole energy subbands calculated 
by means of the SEM for the NWs with thicknesses of 
7 and 15 nm. The subband energies seem to shift with 
temperature because the Fermi level is chosen as an ori- 
gin in the band diagram. Actually, the NW Fermi level 
does shift with temperature. The shift of the Fermi level 
is mainly due to the change in the ratio between the num- 
bers of carriers with opposite polarity with temperature. 
At a low temperature, electrons dominate in the intrinsic 
NWs. For the 15-nm-thick (7-nm-thick) NW, the Fermi 
level overlaps the electron subbands, when the temper- 
ature is less than 190 K (110 K). A number of holes 
increases with temperature; hence, the Fermi level tends 
to the middle of the NW band gap. For the 7-nm-thick 
(15-nm-thick) p— type NWs at the gate voltage of 60 V 
(45 V), the electrons (holes) dominate at a low temper- 
ature; hence, the Fermi level is close to the conduction 
(valence) subband. For the 15-nm-thick p— type NW at 
the excess hole concentration p ex = 5 x 10 18 cm" 3 , the 
Fermi level overlaps the hole subbands when the tem- 
perature is less than 270 K. The separation between the 
energy subbands decreases with the NW thickness. For 



the 7-nm-thick (15-nm-thick) intrinsic NW, the band gap 
decreases from a value of 338 meV (189 meV) at a tem- 
perature of 77 K to a value of 321 meV (180 meV) at a 
temperature of 500 K. For the 7-nm-thick (15-nm-thick) 
p-type NW, the band gap of 493 meV (229 meV) does 
not depend on temperature. The splitting between the 
first and second subbands is much greater than the split- 
ting between the second and third subbands. For the 7- 
nm-thick intrinsic NW, the splitting between the second 
and the third electron subbands increases from a value 
of 14 to 39 meV with temperature, while the correspond- 
ing hole subbands remain very close to each other. For 
the 7-nm-thick (15-nm-thick) p— type NW at the gate 
voltage of 60 V (45 V), the splitting between the sec- 
ond and third electron subbands of 57 meV (18 meV) is 
about 5 (3) times greater than the splitting between the 
corresponding hole subbands of 12 meV (6 meV). This 
effect is due to the difference between the electron and 
hole masses. We recall that, for the nanowires, the elec- 
tron density of states have peculiarities at the edge of 
energy subbands, which leads to the increase in both the 
carrier concentration and the Seebeck coefficient as well 
as the decrease in the thermal conductivity. Hence, for 
the n— type bismuth telluride NWs, the maximal value 
of the figure of merit is achieved at less value of the gate 
voltage in comparison with that for the p— type NWs 
because of the large separation between the second and 
third electron subbands. 

Figure [5] shows that, for the intrinsic NWs with thick- 
nesses of 7 and 5 nm as well as for the p— type NW 
with a thickness of 15 nm at the excess hole concen- 
tration p ex = 5x 10 18 cm" 3 , the NW charge carrier con- 
centration increases with temperature. While, for the 
p— type (n— type) NW with a thickness of 7 nm at 
the excess hole (electron) concentration 5 x 10 18 cm" 3 
and 1 x 10 19 cm" 3 , the NW hole (electron) concentra- 
tion does not depend on temperature. It is equal to 
1.3 x 10 18 cm" 3 (1.4 x 10 18 cm" 3 ) and 2.3 x 10 18 cm" 3 , 
respectively. The temperature dependence of the carrier 
concentration decreases with both an increase in the ex- 
cess hole (electron) concentration and a decrease of the 
NW thickness. Generally, the carrier concentration of 
the considered NWs does not depend on temperature till 
the room temperature. 

Both the SEM and LPE approach are used to obtain 
the data presented in Figs |9"1[T2"1 Figure [5] shows the 
temperature dependence of the Debye screening length 
(Ld) for the 7- and 15-nm-thick NWs. Since the elec- 
tron (hole) gas in the 15-nm-thick intrinsic (p— type) 
NW is degenerate when the temperature is less than 150 
K (see Fig. [Tfa)), the Debye screening length does not 
change with temperature. For the doped NWs, the Debye 
screening length monotonically increases with tempera- 
ture when the Fermi level approaches the middle of the 
band gap. For the 7-nm-thick NWs, the Ljj increases 
by 2.3 times, while it increases by 40 per cent for the 
15-nm-thick p— type NW. In contrast, the temperature 
dependence of the Lu for the intrinsic NWs is not mono- 
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FIG. 7: Temperature dependence of the first three electron and hole subbands (a) for the intrinsic and (b) for the p— type 
bismuth telluride NWs with thicknesses of 7 nm (solid lines) and 15 nm (dashed lines) at the excess hole concentration 
Vex = 5 x 10 18 cm -3 . For the intrinsic (p— type) NWs, the energy spectrum is calculated at gate voltages of 15 and 42 V (60 
and 45 V) for the NW thicknesses 7 and 15 nm, respectively. 
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FIG. 8: (Color online) Temperature dependence of the elec- 
tron and hole concentrations (a) for the intrinsic bismuth tel- 
luride NWs with thicknesses of 7 (solid line), 15 nm (dashed 
line) at gate voltages of 15 V and 42 V, correspondingly, and 
for the p— type bismuth telluride NWs with a thickness of 
15 nm (dashed-dot line) at a gate voltage of 45 V, when 
Pex=5x 10 18 cm" 3 . 



tonic. At a high temperature, the electron concentration 
increases with temperature. Hence, the Debye screening 
length decreases. For the 7- and 15-nm-thick intrinsic 
NWs, the Lb increases by 2 times and by 25 per cent, 
respectively. The Debye screening length also decreases 
with an increase in the excess hole concentration. For 
example, for the 7-nm-thick p— type NW with the excess 
hole concentration p ex = 1 x 10 19 cm -3 , the is less 



than that for the 7-nm-thick NW at p ex — 5 x 10 18 cm -3 
by 25 per cent in the whole temperature range. There- 
fore, the Debye screening length strongly depends on 
temperature for the NWs with a small thickness. 

Figure [TO] shows the temperature dependences of the 
Seebeck coefficient for the NWs with thicknesses of 7 
and 15 nm. For the intrinsic NWs and the 15-nm-thick 
p— type NW, the temperature dependence of the Seebeck 
coefficient (S) is not monotonic, because the electron 
(hole) concentration increases at a high temperature. For 
both the n— and p— type NWs with a thickness of 7 nm, 
at gate voltages of -35 and 60 V, respectively, the Seebeck 
coefficient is negative, because the electrons dominate. 
While, for the 7- and 15-nm-thick p— type NWs with the 
excess hole concentrations of p ex = 1 x 10 19 cm -3 and 
Pex = 5 x 10 18 cm~ 3 , correspondingly, at gate voltages 
of 75 and 45 V, the S is positive owing to the domi- 
nation of holes in the system. For the 7-nm-thick (15- 
nm-thick) intrinsic and 15-nm-thick p— type NWs, the 
maximum absolute values of the Seebeck coefficient 342 
(302) and 267 V/K are achieved at a temperature of 400 
(325) and 390 K, correspondingly. At the room temper- 
ature, for the 7-nm-thick (15-nm-thick) intrinsic NWs, 
the absolute value of the Seebeck coefficient 325 V/K 
(301 V/K) at the gate voltage of 15 V (42 V) is greater 
than a value of 252 V/K (167 V/K) at the zero gate 
voltage by 29 (80) per cent. For the 7-nm-thick p— type 
NW at the excess hole concentration p ex = 5x 10 18 cm~ 3 
(Pex = 1 x 10 19 cm -3 ), the absolute value of S monoton- 
ically increases with temperature by 35 (50) per cent. 
For the 7-nm-thick n— type NW, \S\ is less than that for 
the 7-nm-thick p— type NW by 5 per cent in the whole 
temperature interval. The monotonic behavior of the 
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FIG. 9: (Color online) Temperature dependence of the Debye screening length (a) for the intrinsic bismuth telluride NWs 
with thicknesses of 7 (solid line) and 15 nm (dashed line) at gate voltages of 15 V and 42 V, correspondingly, and (b) for the 
p— type bismuth telluride NWs with thicknesses of 7 nm (solid line) and 15 nm (dashed line) at gate voltages of 60 and 45 V, 
respectively, when p ex = 5 x 10 18 cm -3 . At V g = —35 V (V g = 75 V), the dashed-dotted line (dotted line) corresponds to the 
Debye screening length for the n— type (p— type) bismuth telluride NW with a thickness of 7 nm when n ex = 5x 10 18 cm~ 3 
(p ex = 1 x 10 19 cm" 3 ). 




FIG. 10: (Color online) Temperature dependence (a) of the Seebeck coefficient (S) for the intrinsic bismuth telluride NWs 
with thicknesses of 7 nm (solid line) and 15 nm (dashed line) at gate voltages of 15 and 42 V, correspondingly, and (b) of the 
absolute values of S for the p— type bismuth telluride NWs with thicknesses of 7 nm (solid line) and 15 nm (dashed line) at gate 
voltages of 60 and 45 V, respectively, when p ex = 5 x 10 18 cm -3 . At Vg =-35 V (Vg =75 V), the dashed-dotted line (dotted 
line) corresponds to the S for the n— type (p— type) bismuth telluride NW with a thickness of 7 nm when n ex = 5 x 10 18 cm -3 
{ Pex = 1 x 10 19 cm- 3 ). 



Seebeck coefficient as a function of temperature results 
from the weak local variation of the electron concentra- 
tion, when the Debye screening length is greater than 
the NW thickness (see Fig. [5Jb)). At the room tempera- 
ture, for the 7-nm-thick p— type (n— type) NWs with the 



excess hole (electron) concentration p ex = 5 x 10 18 cm~ 3 
(n ex = 5x 10 18 cm~ 3 ) and p ex — 1 x 10 19 cm" 3 , the ab- 
solute value of the Seebeck coefficient 319 (302) and 276 
V/K at gate voltage of 60 (-35) and 75 V, respectively, is 
greater than a value of 206 (187) and 142 V/K at the zero 
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gate voltage by 54 (61) and 95 per cent, correspondingly. 
For the 15-nm-thick p— type NW at the room tempera- 
ture, the Seebeck coefficient of 196 V/K increases by 25 
per cent when applying the gate voltage of 45 V. Hence, 
the dependence of the Seebeck coefficient on tempera- 
ture is diminished with a decrease in the NW thickness, 
while the absolute value of S increases under the applied 
electric field. 

Figure fTTI presents the temperature dependence of the 
thermal conductivity of the bismuth telluride NWs with 
thicknesses of 7 and 15 nm. For comparison, we depict 
the temperature dependence of the thermal conductivity 
for the NWs with the same thicknesses at a zero gate 
voltage. For the 7- and 15-nm-thick intrinsic NWs, the 
thermal conductivity at the gate voltages of 15 and 42 
V, respectively, is less than that at the zero gate voltage 
when the temperature is greater than 310 and 240 K. At 
the room temperature, for the 15-nm-thick (7-nm-thick) 
intrinsic NW, the thermal conductivity of 0.434 W/Km 
(0.274 W/Km) at the gate voltage of 40 V (15 V) is less 
(greater) compared to the value of 0.727 W/Km (0.254 
W/Km) at the zero gate voltage by 40 (8) per cent. The 
temperature dependence of the thermal conductivity is 
not monotonic because its electron (hole) part increases 
with temperature, while its lattice part decreases. At 
a high temperature, for the intrinsic NWs, the temper- 
ature dependence of the thermal conductivity is much 
stronger in comparison with that for the doped NWs. 
The increase in the thermal conductivity of the intrinsic 
NWs arises from the additional thermal transport chan- 
nel produced by the electron-hole pairs (see Eq. [101) • The 
electric field bias enables the charge carrier concentra- 
tion to increase. As a result, for the 7-nm-thick (15-nm- 
thick) intrinsic NW at the gate voltage of 15 V (42 V), 
the minimum value of the thermal conductivity of 0.268 
W/Km (0.352 W/Km) at T = 275 K (247 K) is greater 
than a value of 0.195 W/Km (0.259 W/Km) at the zero 
gate voltage by 36 per cent, when T = 247 K (190 K). 
For the doped NWs, the applied lateral electric field de- 
creases the thermal conductivity in the whole temper- 
ature range. For the 7-nm-thick (15-nm-thick) p— type 
NW at the gate voltage of 60 V (45 V), the minimum 
value of the thermal conductivity of 0.219 W/Km (0.277 
W/Km) at T = 360 K (275 K) is less than the value of 
0.307 W/Km (0.384 W/Km) at the zero gate voltage by 
27 per cent when T = 246 K. At the room temperature, 
for the 7-nm-thick (15-nm-thick) p— type NW at the ex- 
cess hole concentration p ex = 5 x 10 18 cm -3 , the thermal 
conductivity of 0.222 W/Km (0.278 W/Km) at the gate 
voltage of 60 V (45 V) is less than a value of 0.315 W/Km 
(0.398 W/Km) at the zero gate voltage by 30 per cent. 
For the 7-nm-thick p— type NW, the temperature depen- 
dence of the thermal conductivity is identical to that for 
the n— type NW. Therefore, the applied electric field per- 
mits decreasing the NW thermal conductivity. 

Figure [T^] shows the temperature dependence of the 
figure of the merit for the intrinsic and doped NWs with 
thicknesses of 7 and 15 nm. At the gate voltage of 15 



V (42 V), the intrinsic 7-nm-thick (15-nm-thick) NW ex- 
hibits the maximum value of the figure of merit as high 
as 3.8 (2.4) at the temperature of T = 360 K (300 K). As 
a result of the great increase in the thermal conductivity 
and decrease in the Seebeck coefficient for the intrinsic 
7-nm-thick (15-nm-thick) NW within the high temper- 
ature range, the figure of merit decreases rapidly with 
temperature when T > 360 K (300 K). For the 7-nm- 
thick (15-nm-thick) p— type NWs with the excess hole 
concentration p ex = 5 x 10 18 cm -3 , the maximum value 
of the figure of merit is 5.12 (1.83) at a temperature of 
500 K (360 K), when the gate voltage is 60 V (45 V). As a 
function of temperature, the figure of merit for the 7-nm- 
thick n— type NW resembles that for the p— type NW. 
Their figure of merit is about twice as large as that at the 
zero gate voltage. For the 7-nm-thick doped NWs, the 
monotonic dependence of the figure of merit on temper- 
ature is attributed to the slight temperature dependence 
of the Seebeck coefficient and the thermal conductivity. 
As a result of decrease of the Seebeck coefficient with an 
increase in the excess hole concentration, the figure of 
merit decreases. For example, at the temperature of 500 
K, for the 7-nm-thick p— type NW with the excess hole 
concentration p e2 ; = 1 x 10 19 cm" 3 , the figure of merit of 
3.61 is less than that corresponding to the NW with the 
excess hole concentration n ex = 5 x 10 18 cm" 3 by 30 per 
cent. For the intrinsic NWs, the extreme point of the 
figure of merit is close to the room temperature. Both 
the confinement effect and the excess holes (electrons) 
produce a shift of the extreme point toward a high tem- 
perature. Hence, for the intrinsic NWs, the electric field 
effect gives rise to an increase in the figure of merit by an 
order of magnitude; while, for the doped NWs, the figure 
of merit can increase as much as twice. 



VI. CONCLUSIONS 

In this paper, we showed that a combination of the 
quantum confinement effect in the bismuth telluride 
quantum wires with a properly applied electric bias can 
increase the thermoelectric power conversion by an order 
of magnitude. For the considered bismuth telluride NWs, 
the electron-electron exchange-correlation interaction is 
negligible because of the small electron (hole) concen- 
tration. At room temperature, the exchange-correlation 
energy does not exceed the electron thermal energy. The 
EFE can effectively control the type of the NW conduc- 
tivity by adjusting the charge carrier concentration. It 
compares favorably with the doping effect because the 
number of scattering centers is not changed under an the 
electric field bias. The variation in the charge carrier 
concentration leads to a self-consistent change of the size 
quantization effect. For bismuth telluride NWs, the elec- 
tron effective mass is less than the hole effective mass. 
As a result, the splitting between the electron subbands 
is greater than that between the hole subbands. In ad- 
dition, for bismuth telluride material, the electron mo- 
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FIG. 11: (Color online) Temperature dependence of the thermal conductivity (a) for the intrinsic bismuth telluride NWs with 
thicknesses of 7 nm (solid line) and 15 nm (dashed line) at gate voltages of 15 and 42 V, correspondingly, and (b) for the 
p— type bismuth telluride NWs with thicknesses of 7 nm (solid line) and 15 nm (dashed line) at the gate voltages of 60 and 45 
V, respectively, when p ex = 5 x 10 18 cm -3 . At V g = —35 V (V g = 75 V), the dashed-dotted line (dotted line) corresponds to 
the thermal conductivity for the n— type (p— type) bismuth telluride NW with a thickness of 7 nm when n ex = 5 x 10 18 cm -3 
{pex = 1 X 10 19 cm -3 ). At the zero gate voltage, the dependence of the thermal conductivity is shown by the dashed-dotted 
(filled triangles) and dotted lines (filled squares) for the intrinsic (p— type) NWs with thicknesses of 7 and 15 nm, respectively, 
in Fig. Illf a) fFig. Illf b)). The filled squares (triangles) correspond to the data obtained in the LPE approach for the NW with 
a thickness of 15 nm (7 nm) (see Fig. Illf a)). 




FIG. 12: (Color online) Temperature dependence of the figure of merit (a) for the intrinsic bismuth telluride NWs with 
thicknesses of 7 nm (solid line) and 15 nm (dashed line) at gate voltages of 15 and 42 V, correspondingly, and (b) for the 
p— type bismuth telluride NWs with thicknesses of 7 nm (solid line) and 15 nm (dashed line) at the gate voltages of 60 V and 
45 V, respectively, when p ex = 5 x 10 18 cm -3 . At V g = —35 V (V g = 75 V), the dashed-dotted line (dotted line) corresponds 
to the figure of merit for the n— type (p— type) bismuth telluride NW with a thickness of 7 nm when n ex = 5 x 10 18 cm -3 
{Pex = 1 X 10 19 cm -3 ). At the zero gate voltage, the dependence of the figure of merit is shown by the dashed-dotted (filled 
triangles) and dotted lines (filled squares) for the intrinsic (p— type) NWs with thicknesses of 7 and 15 nm, respectively, in 
Fig. I12f a) (Fie. 112T b)). The filled squares (triangles) correspond to the data obtained in the LPE approach for the NW with a 
thickness of 15 nm (7 nm) (see Fig. 112( a)). 



16 



bility is about twice as high as the hole mobility. It in- 
dicates that the thermoelectric properties of the n— type 
NWs are superior to those of the p— type NWs. The 
dependence of the Seebeck coefficient on temperature is 
diminished with a decrease in the NW thickness, while 
the absolute value of S increases under the applied elec- 
tric filed. In dependence on the NW thickness and ex- 
cess hole (electron) concentration, the electric field bias 
enables the Seebeck coefficient to increase by 30-90 per 
cent. For the 7-nm-thick (15-nm-thick) intrinsic and 15- 
nm-thick p— type NWs, under the applied electric field, 
the maximum absolute values of the Seebeck coefficient 
of 342 (302) and 267 V/K are achieved at a tempera- 
ture of 400 (325) and 390 K, correspondingly. The con- 
finement effect on the electron density of states leads to 
an increase in the NW Seebeck coefficient, while the in- 
crease in the carrier concentration causes the decrease in 
S. The applied electric field leads to a decrease in the 
electronic contribution to the NW thermal conductivity. 
For example, at the room temperature, for the 15-nm- 
thick p— type (intrinsic) NW with the excess hole con- 
centration p e 



5 x 10 18 cm 3 the thermal conductiv- 



ity of 0.278 W/Km (0.426 W/Km) at gate voltage of 45 
V (42 V) is less compared to the respective value of 0.398 
W/Km (0.727 W/Km) at the zero gate voltage by 30 (40) 
per cent. For the 7-nm-thick (15-nm-thick) p— type NWs 
with the excess hole concentration p ex = 5 x 10 18 cm~ 3 , 
the maximum value of the figure of merit is 5.12 (1.83) at 
a temperature of 500 K (360 K), when the gate voltage is 
60 V (45 V). At the room temperature, for the 7-nm-thick 
(15-nm-thick) p— type NW with the excess hole concen- 
tration of p ex = 5 x 10 18 cm -3 , the maximum value of 
the figure of merit of 3.1 (1.5) is greater than the corre- 
sponding value ZT = 1.65 (1.19) at the zero gate voltage 
by 90 (25) per cent. For the 7-nm-thick (15-nm-thick) in- 
trinsic NW, the maximum value of the figure of merit is 
as high as 3.4 (2.4) at the gate voltage of 15 V (42 V). It 
is an order of magnitude greater than the corresponding 
value of the figure of merit at the zero gate voltage. 



local variation of the electron and hole concentrations of 
each valley in a nanowire a o 22 i 34 



n 1D (x, z) 



N, 



ID n 



a x a z 



E E $ -i/2 



2 = 1 



± 



eip(x, z) 



k B T 



(Al) 

Here, the factor = {2m e y h k B T /nh 2 ) 1/2 denotes the 



effective density of electron (hole) states. Sign "+" ("-") 
corresponds to electrons (holes). The reduced chemical 

potentials rjffi for an electron subband and a hole sub- 
band are defined as 
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For a non-degenerate semiconductor material, if the 
applied electric filed is small enough compared to 
the difference between the Fermi energy and the bot- 
tom (top) of the conduction (valence) band such 



Ef — E r 



E 



n . / 



>> \etp(x, z)\, then the electron 

(hole) concentration n\D (pid) can be expanded in Tay- 
lor series in terms of eLp/ksT . Using this expansion, the 
right hand side of Eq. ([2]) can be linearized with regard to 
the electrostatic potential <p. Hence, Poisson's equation 
can be rewritten in the linearized form as 



APPENDIX A: LINEARIZATION OF POISSON'S 
EQUATION 
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In the semi-classical approximation, to take into ac- 
count the electron-screening effects, one can write the 



Here, the Dcbye screening length is defined by 
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Since the electric field is directed along the z axis and of the intrinsic NWs, the variables in Eqs. |T]) and (|A4j) 
the exchange-correlation energy V xc is negligible in case can be separated. Therefore, the analytical solution of 
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the boundary problem (J6j> , |[7J), and (|A4[) reads 
L D e S i0 2 smh{z / L D ) 



<p{z) = V g 



dsi0 2 CBi 2 Te 3 cosh(a z /L D ) 



(A6) 



The solution of the Schrodinger equation for electrons 
reads 



& n j(x,z) = \ — sm(nlx/a x )Z n (z). (A7) 
a x 

The unknown function Z n (z) satisfies the following equa- 
tion 



h 2 d 2 Z n 
2m z dz 2 



eip{z)Z n (z) = X n Z n (z), 



(A8) 



and the corresponding boundary conditions 
Z n (0) = Z n {a z ) = 0. The electron subband energy 
is given by the expression l = A„ + h 2 ir 2 l 2 /2m x a 2 . 
So, we recast the self-consistent system of differential 
equations (Q} and ([2]) into the simplified self-consistent 
system of one differential equation (|A8|) and two mathe- 
matical expressions (|A5|) and (|A6j) taking into account 
charge neutrality equation Q. To solve the last system 
of equations, we apply the following algorithm: STEP 1: 
Initial guess of the Debye screening length STEP 
2: Calculation of the electrostatic potential ip( k '{z) 

by means of Eq. JA6]) using L$ . STEP 3: Solving 

(fc) 



Eq. jA8]) for the eigenvalues A„ using tp {k >(z). STEP 
4: Solving Eq. ^ for the Fermi energy E 



(k) . „(fc) 

using E K n j. 



STEP 5: Calculation of L 



(fc+i) 



by means of Eq. 



using E K '. 



STEP 6: IF 



-(fe+i) 

J D 



L 



/L 



(fc+i) 

D 



> 5 



THEN L 



_ . - [L ( ^ +1) + LP) /2 AND GOTO STEP 2 

ELSE STOP. The solution of Eq. {AS} can be obtained 
analytically and numerically. To find the analytical 
solution, one transforms Eq. (|A8[) into the modified 
Mathieu equation i 23 ' 24 Therefore, in case of intrinsic 
semiconductor NWs, the linearized Poisson equation 
(LPE) approach allows us to simplify solving the self- 
consistent system of Schrodinger and Poisson equations 
significantly. To solve equation (|A8|) numerically, we can 
use the spectral method or spectral element method . 19 i 35 



APPENDIX B: SPECTRAL METHOD WITH 
TRIGONOMETRIC BASIS FUNCTIONS 



The function Z(z) in Eq. (|A8[) can be decomposed in 
terms of the trigonometric basis functions 



JV 
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fc=i 



Here, we set the reduced variable t = z/Ljj which ranges 
in the interval [0, r\. The upper limit of integration is 
T = a z /Lo- The basis consists of the wave func- 

tions $fc(i) = \/2/r sin(7rfci/r) of an electron in the one- 
dimensional Infinite Square well potential. To find the 
unknown constants Ck we substitute expression (|B1() into 
Eq. (|A8[) . multiply both sides of the equation by the con- 
jugated function (i) and integrate from to r. Finally, 
we obtain a system of linear algebraic equations for N 
unknowns Cfc. The eigenvalues A„ can be calculated by 
means of the secular (characteristic) equation 
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where m and n — 1,2, ... N and A = 2\m z L 2 D /h 2 . We 
take N — 30 here. The integral in Eq. (]F32|) is calculated 
analytically when the electrostatic potential ip(t) is given 
by expression (|A6[) . In the general case, if the potential in 
Schrodinger equation |T]) presents a complex function, we 
can solve this equation by using the more sophisticated 
method, that is, spectral element method described in 
the next subsection. 



APPENDIX C: SPECTRAL ELEMENT METHOD 
BASED ON THE GAUSS-LOBATTO-LEGENDRE 
POLYNOMIALS 

In the SEM approximation, the component of the elec- 
tron wave function and the electrostatic potential are ex- 
panded in terms of the field variables 19 
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where the reduced variable z ranges in the interval 
— 1 < z < 1. The field variables in the Schrodinger 
and Poisson's equations are represented by the Gauss- 
Lobatto-Lelendre polynomials bj(z) of the Nth order 



bJz) 



(l-z 2 )L' N (z) 



NiN+^LvizjXzj-zy 



(C3) 



where Ljv(z) is the Nth order Legendre polynomial. 
The jth zero Zj is a root of the equation L' N (zj) = 0. 



We suppose 7V=15 here. To obtain a weak form of 
the Schrodinger and Poisson's equations, we substitute 
expressions ([CTj) and JC2J into the modified Eqs. (1]) 
and @, multiply both sides of these equations by 
the function bj(z) and integrate from -1 to 1. Fur- 
thermore, we apply the Gaussin-Legendre quadrature 
wi = 2/(1 — zf)L 2 N+1 (zi) to evaluate the integration and 
finally obtain the discrete Schrodinger and Poisson's 
equations as 



f 1 dbjjzi) dbjjzi) 2Jm z a 2 \ 
Z^l^ Wl \ ~J~ Hi dT~ — W — \- V ^ Zl > + E c- e( P( z i)\ bi(zi)b 3 (zi) <> Z{Zj) 
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1=0 



where i = 0, 1, ...,N , X = 2Xm z a 2 /h 2 , and J is the 
Jacobian. Eq. (|C4|) can be solved by a generalized 
eigenvalue solver. However, because the mass matrix 
in Eq. (|C4[) is diagonal, this equation can be converted 
into a more efficient regular eigenvalue problem.— 
Eq. (|C5|) can be solved by the Newton-Raphson method 
using the predictor-corrector approach i 25 ' 26 i 27 To solve 
the self-consistent system of Schrodiger and Poison's 
equations (|C4j) and (|C5[) . we use the following algorithm: 
STEP 1: Suppose the electron-electron exchange- 
correlation energy is Vic\z)=0. STEP 2: Solving 
Poisson's equation @ with boundary conditions (O 
and ([7]) for the electrostatic potential ip(°'(z), when 
p = 0. STEP 3: Solving Schrodinger equation (|C4|) for 

the eigenvectors Z^\zj) and eigenvalues taking 
into account the electrostatic potential ip( k \zj) and the 
electron-electron exchange-correlation energy Vic } (,?,). 
STEP 4: Solving Eq. (g]) for the Fermi energy using 



E^j. STEP 5: Solving Eq. ([C5} for <^ (fe) (%) by means 
of the Newton-Raphson metod using the predictor- 
corrector approach. STEP 6: Calculation of the 
electron-electron exchange-correlation energy Vxc\zj) 
by means of Eq. ([3]) using the Fermi energy E^ , energy 
spectrum E^J , and wave functions 4^) (zj). STEP 7: IF 

Ef= - M fe ^)l]VEfiol^ fc) (%)l 2 > s 

THEN (p( k+ V( Zj ) := ip {k) {zj) AND GOTO STEP 3 
ELSE STOP. 
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